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ABSTRACT 


i 

A  coplanar  Schottky  variable  phase  shifter  is  studied.  This 
device  has  a  simple  configuration  and  Is  suitable  for  monolithic  inte¬ 
gration.  The  device  uses  the  electronic  variability  of  the  depletion 

layer  thickness  of  the  Schottky  contact  under  the  coplanar  waveguide  to 

» 

change  the  amount  of  the  phase  shift.  A  periodic  structure  Is  Introduced 
to  reduce  a  loss  caused  by  the  existing  semiconductor  material  which  has 
a  finite  conductivity.  An  Integral  equation  Is  formulated  and  solved  by 
the  point-matching  method.  The  analysis  allows  one  to  predict  the  phase 
shift  and  the  attenuation  In  the  device.  It  Is  shown  that  such  a  device 
can  operate  at  millimeter-wave  frequencies  with  a  very  small  loss.  A 
simple  measurement  confirmed  these  predictions. 
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CHAPTER  1  :  INTRODUCTION 


Transmission  lines  constructed  on  semiconductor  substrate  were 
first  studied  to  investigate  the  behavior  of  the  signal  travelling  in  a 
semiconductor  chip  in  a  high  speed  computer.  Several  analyses  based  on 
parallel-plate  models  have  been  published,  and  it  has  been  shown  that 
such  structures  exhibit  slow-wave  characteristic  under  certain  circum¬ 
stances  [1-3].  The  electrical  length  of  such  a  structure  could  be  forty 
to  fifty  times  longer  than  the  physical  length.  This  means  that  the 
effective  dielectric  constant  could  become  one  thousand  or  more.  The 
waveguide  structure  usually  contained  two  different  media:  an  insulator 
and  a  doped  semiconductor,  which  is  often  refered  to  as  an  MIS 
(  metal-insu  ator-semiconductor  )  structure.  The  slow-wave  character¬ 
istic  is  caused  by  the  phenomenon  that  the  electric  and  magnetic  energies 
are  stored  in  different  regions  of  the  waveguide.  The  slow-wave  phenome¬ 
non  is  interesting  In  microwave  and  millimeter-wave  (  mm-wave  )  applica¬ 
tions  such  as  delay  lines  and  phase  shifters.  These  transmission  lines 
have  either  MIS  configurations  or  Schottky-contacts  on  semiconductor 
substrates.  In  the  latter  case,  the  depletion  layer  formed  under  the 
Schottky-contact  metal  acts  like  an  insulator  region  since  it  has  a  small 
conductivity  compared  to  the  surrounding  doped  semiconductor  region. 
The  propagation  speed  depends  greatly  on  the  thickness  of  the  insulated 
layer  under  the  metal.  Therefore,  the  bias  dependence  of  the  Schottk- 
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y-contact  depletion  layer  thickness  can  be  used  advantageously  to  con¬ 
struct  an  electronically  variable  phase  shifter  [4-6]. 

In  practice,  various  planar  structures  can  be  used  to  realize 
the  slow-wave  structure,  such  as  microstrip  line  and  coplanar  waveguide. 
Between  these  structures,  a  coplanar  waveguide  seems  to  be  preferable 
for  various  applications  because  other  devices  can  be  connected  In 
series  or  parallel.  For  this  reason,  this  discussion  Is  limited  to  the 
MIS  and  Schottky  coplanar  waveguide.  These  structures  are  also  suitable 
for  applications  In  monolithic  microwave  and  millimeter-wave  integrated 
circuits. 

A  large  problem  of  these  devices  is,  however,  their  Inherently 
large  attenuation  caused  by  the  finite  conductivity  of  the  semiconductor 
substrate.  To  obtain  slow-wave  propagatin  up  to  a  high  frequency  range, 
It  Is  necessary  to  raise  the  conductivity  of  the  semiconductor  substrate 
to  a  large  value,  typically  10*’  to  10s  S/m.,  which  causes  high  atten¬ 
uation  at  high  frequencies.  Not  many  investigators  have  tried  to  solve 
this  problem.  There  were  a  few  papers  dealing  with  the  reduction  of  the 
loss  by  replacing  the  semiconductor  material  by  periodic  metal  strips 
[7-8].  The  slow-wave  propagation  occurlng  In  the  structures  Is  caused  by 
the  periodicity.  Therefore,  this  structure  only  solves  the  attenuation 
problem  for  a  fixed  phase-delay  line,  and  can  not  be  applied  to  the 
Schottky-contact  line  with  the  property  of  the  electronic  variability  of 
the  depletion  layer  thickness. 

In  Part  I  of  this  dissertation,  a  new  structure  is  presented  to 
solve  this  problem:  an  MIS  and  Schottky  periodic  coplanar  waveguide.  As 
Fig.  1  shows,  the  proposed  structure  Introduces  periodic  lossless  sec- 


conductors 


slots 


Figure  1  Schematic  view  of  the  coplanar  waveguide 
constructed  on  the  periodically  doped 
semiconductor  substrate 
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tlons  Into  the  MIS  or  Schottky  coplanar  waveguide.  In  this  way,  the 
total  attenuation  of  the  device  can  be  reduced  since  it  has  fewer  lossy 
semiconductor  regions.  Also,  the  periodicity  of  the  structure  provides 
slow-wave  propagation,  and  the  total  phase  delay  value  can  be  improved. 

Part  I  covers  both  theoretical  and  experimental  study  of  such  a 
structure.  The  theoretical  analysis  Includes  the  MIS  coplanar  waveguide 
using  an  Integral  equation  method  and  also  the  analysis  of  the  periodic 
structure  by  conventional  transmission  line  theory.  The  experiments 
were  conducted  performed  to  confirm  the  theory.  Several  models  of  the 
proposed  structure  were  constructed,  and  tested  at  VHF  and  UHF  frequen¬ 
cies.  The  results  confirm  the  applicability  of  the  theory.  In  the  last 
chapter  of  Part  I,  optimum  conditions  for  the  operation  of  the  coplanar 
Schottky  variable  phase  shifter  are  discussed  using  the  theoretical 
results.  It  will  be  shown  that  such  a  device  can  be  used  at  mm-wave  fre¬ 
quencies  without  the  disadvantage  of  significant  attenuation. 


CHAPTER  2  :  ANALYSIS  OF  MIS  COPLANAR  WAVEGUIDE 


In  this  chapter,  the  propagation  characteristics  of  the  MIS 
coplanar  waveguide  are  discussed.  An  elaborate  method  of  analyzing  this 
structure  is  developed  which  is  capable  of  predicting  the  frequency 
dependence  of  the  propagation  characteristics.  From  the  result  of  this 
analysis,  the  frequency  region  of  the  slow-wave  propagation,  the  value 
of  the  slow-wave  factor,  the  attenuation  constant,  and  the  character¬ 
istic  impedance  can  be  determined.  Since  the  MIS  coplanar  waveguide  Is  a 
special  case  of  the  MIS  periodic  coplanar  waveguide,  the  results  of  this 
analytical  method  are  used  in  the  following  chapters  to  determine  the 
characteristics  of  the  periodic  structure  as  well. 


METHOD  OF  ANALYSIS 

Several  previous  studies  have  shown  the  applicability  of  diverse 
techniques  to  the  analysis  of  MIS  coplanar  waveguide.  The  spec¬ 
tral-domain  method,  mode-matching  method,  and  finite  element  method  have 
all  been  used  [9-15].  The  convergence  of  the  solutions  in  these  methods 
Is  dependent  on  the  actual  field  distribution  in  the  MIS  coplanar  wave¬ 
guide  and  the  choice  of  basis  functions  in  the  analysis.  In  the  MIS 
coplanar  waveguide,  the  field  around  the  slot  region  is  deformed  by  the 
highly  conductive  doped  region  located  near  the  coductors.  Therefore, 
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the  spatial  gradients  of  the  field  components  are  large  near  the  edges  of 
the  conductors,  and  hence,  many  basis  functions  are  needed  In  that  region 
to  express  them  [15]. 

Yamashlta  Introduced  a  hybrid-mode  analysis  technique  to  solve 
lossless  planar  structures  with  various  configurations  [16].  In  his 
technique,  an  Integral  equation  Is  first  constructed  with  respect  to 
unknown  field  components,  and  Is  then  solved  by  non-uniform  discretiza¬ 
tion.  Matching  points  are  taken  so  that  more  points  are  located  near  the 
edges  of  the  conductors  to  make  the  convergence  of  the  solutions  faster. 
This  technique  Is  especially  well  suited  for  the  structure  under  consid¬ 
eration,  because  larger  spatial  gradients  of  the  field  components  near 
the  edges  of  the  conductors  are  expected  compared  to  the  conventional 
lossless  structures.  Another  advantage  of  this  technique  Is  that  It  can 
easily  be  used  to  calculate  the  characteristic  Impedance  of  the  struc¬ 
ture  by  choosing  appropriate  field  components  as  the  unknown  functions 
of  the  integral  equation.  Also,  this  method  Is  a  full-wave  analysis,  and 
Is  used  to  obtain  the  frequency  dependence  of  the  characterises.  This  Is 
Important  since  the  slow-wave  phenomenon  of  the  MIS  structure  strongly 
depends  on  frequency. 

For  convenience,  the  cross  section  of  the  MIS  coplanar  waveguide 
Is  modeled  as  In  Fig.  2.  The  substrate  Is  divided  into  three  layers  each 
of  uniform  thickness:  the  Insulated  layer,  the  doped  layer,  and  the 
semi-insulated  layer  (  Regions  2,  3,  and  4  In  the  figure, 
respectively  ).  These  regions  are  assumed  to  be  linear,  homogeneous, 
and  Isotropic.  That  Is,  material  properties  are  uniform  In  each  region, 
and  are  characterized  by  a  single  dielectric  constant  and  conductivity. 


Figure  2  Cross-sectional  view  of  the  analytical 


model  of  the  MIS  coplanar  waveguide 


The  conductivity  of  the  material  is  Included  in  the  imaginary  part  of  the 
dielectric  constant.  The  metal  layer  is  assumed  to  be  infinitesimally 
thin  and  to  be  perfectly  conductive.  Because  the  structure  is  symmetric 
and  we  are  interested  only  in  the  dominant  mode,  which  has  even  symmetry, 
we  can  place  a  vertical  magnetic  wall  at  the  center  (  x=0  )  and  consider 
only  the  right  half  of  the  structure.  The  tangential  components  of  the 
magnetic  field  vanish  at  the  magnetic  wall.  An  electric  wall  (  perfect 
conductor  )  is  then  placed  at  the  far  right  (  x=w,  w»b  )  for  analytical 
convenience.  This  way,  we  can  generate  a  discrete  Fourier  series  sol¬ 
ution. 

The  analysis  is  based  on  the  electric  and  magnetic  vector  poten- 
tlals,  F  and  A,  respectively  [17].  These  potentials  satisfy  the  follow¬ 
ing  equations: 


7  x  v  x  A 


1 

7"V*a 

J  WE 


(  1  ) 


V  X  V  X  F  - 


Own 


(  2  ) 


where  k  is  the  complex  wave  number,  w  is  angular  frequency,  e  is  the  com¬ 
plex  dielectric  permittivity,  v  is  the  magnetic  permeability,  and  #  and 

Q 

are  arbitrary  scalar  functions.  In  our  case,  it  is  convenient  to 


choose 
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A  =  1>  y 


(  3  ) 


F  =  #  y 


(  4  ) 


where  <p  and  <f  are  the  scalar  functions  which  represent  the  wave  trans¬ 
verse  electric  (  TE  )  wave  and  transverse  magnetic  wave  (  TM  )  to  the 
y-direction,  respectively,  and  y  Is  a  unit  vector  in  the  y-direction. 
The  electric  and  magnetic  field  vectors  can  be  derived  from  these  quanti¬ 
ties: 


->  -►  1 

E  =  -  V  x  F  +  -  Vx?xA 

jwe 


(  5  ) 


-*■  -*■  1  -*■ 

H  =  V  x  A  +  -  V  x  v  x  F 

jwv 
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where  E  and  H  are  the  electric  and  magnetic  field  vectors,  respectively. 
These  potentials  are  expanded  in  terms  of  their  eigenfunctions,  which 
should  be  either  sinusoidal  or  exponential  functions  and  must  satisfy 
the  boundary  conditions  on  the  vertical  magnetic  and  electric  walls.  The 
potentials  in  each  region  can  be  expressed  as  follows: 


(  Region  1  )  (  7  ) 


I  cosBnx  (  C„  sin«2ny  +  C'„  coS«2ny  ) 
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(  Region  2  )  (  8  ) 
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(  Region  3  )  (  9  ) 

J  c°spnx  [  Gn  s1na4n(y+d1+d2)  +  G*n  cosa^Cy+dj+d^  ] 
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(  Region  4  )  (  10  ) 
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(  Region  5  )  (  11  ) 


The  z-dependence  Is  omitted  from  each  of  the  above  expressions  and 
assumed  to  be  of  the  form: 


-j*z 
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where  Z  is  the  propagation  constant  in  the  z-dlrectlon  and: 
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The  summations  over  n  range  from  one  to  infinity,  and  An,  Bn>  etc.  are 
constants  to  be  calculated  from  the  boundary  conditions.  The  potentials 
In  Regions  1  and  5  are  chosen  so  that  the  field  decays  exponentially  as  y 
approaches  ±«  In  each  region. 

The  coefficients  of  the  potentials  are  adjusted  to  satisfy  the 
boundary  conditions  at  y=-dj,  -(d^+dg),  and  -(dj+d2+d3).  These  are  eas¬ 
ily  adjusted  by  comparing  the  coefficients  of  the  terms  of  the  same  value 


of  n  In  the  summation  In  separate  expressions,  since  the  associate  eigen¬ 
functions  have  the  same  variation  In  the  x-dlrection.  Doing  this  elimi¬ 
nates  a  number  of  the  unknown  coefficients,  and  the  potentials  in  the 
Insulated  region  (  Region  2  )  can  then  be  written  as 


e2°ln 

*2  ~  -  1  -  An  cos6nx  <  51  "W  +  Pcn  cos*2ny  > 

"  ‘la2n 

*2  *  1  B„  s,n6»x  <  Pdn  s,"“2ny  *  cos“2ny  > 
n 

(  12  ) 

Pcn  and  Pcjn  are  the  constants  determined  by  this  process.  Derivation  of 
these  values  Is  explained  In  Appendix  A. 

The  x  and  z  components  of  the  electric  field  are  required  to  be 
continuous  at  the  plane  where  conductors  are  placed  (  y=0  )  and  to  be 
zero  on  the  conductors  (  y=0  and  0<x<a,  b<x<w  ),  and  the  x  and  z  compo¬ 
nents  of  the  magnetic  field  to  be  continuous  in  the  slot  region  (  y=0  and 
a<x<b  ).  This  process  Is  summarized  as  follows: 


At  y=0, 
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After  satisfying  these  conditions,  standard  Fourier  analysis  is  used  to 
derive  the  Integral  equation: 

2  [  Pn  J  f(x')  s1n0nx'  dx'  +  Qn  /  g(x')  cos^x'  dx'  ]  cos3nx  =  0 

ns  c 

0<x<a,  b<x<w 

T  [  Rn  T  f(x')  sin3nx'  dx*  +  Pp  /  g(x')  cos&nx'  dx*  ]  sin&nx  =  0 
ns  c 

a<x<b 
(  H  ) 

where  labels  s  and  c  represent  the  Integration  range  over  the  slot  and 
conductors,  respectively.  f(x)  is  proportional  to  the  x  component  of  the 
electric  field  in  the  slot  region  and  g(x)  is  proportional  to  the  current 
density  on  the  conductors: 


a<x<b 


9(x)  -  Hxl  -  Hx2  0<x<a,  b<x<w 


(  15  ) 


The  expressions  of  the  coefficients  P.,,  Q  .  and  R  are  also  given  in 

n  n  n  3 

Appendix  A.  The  main  reason  to  choose  and  Jz  as  unknown  functions  in 
the  Integral  equations  Is  that  It  will  be  easy  to  calculate  the  charac¬ 
teristic  impedance  of  the  structure  after  we  obtain  the  solution. 
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The  integral  equation  (  14  )  is  solved  by  the  point  matching 
method.  Matching  points  are  selected  non-uniformly  along  the  y=0  plane 
to  make  the  convergence  of  the  solution  faster.  More  points  are  selected 
near  the  edges  of  the  conductors  to  accommodate  the  gradients  in  the 
field  components.  The  unknown  functions  Ex  and  J2  are  assumed  to  be  con¬ 
stant  between  two  adjacent  matching  points.  Then  a  system  of  linear 
equations  can  be  obtained.  This  is  written  in  matrix  form: 

Z  u  =  0 

(  16  ) 

where  Z  is  a  two-dimensional  matrix  with  the  size  of  N  xN  ,  where  N  is 

P  P  P 

the  number  of  the  matching  points,  and  a  column  vector  u  contains  unknown 
function  values  f(x)  and  g(x)  between  matching  points.  A  nontrivial  sol¬ 
ution  Is  obtained  only  If  the  matrix  Z  is  singular.  Therefore,  we  seek 
this  condition  by  varying  the  value  of  the  propagation  constant.  This  is 
a  complex  value  since  the  equation  is  complex.  Both  phase  and  atten¬ 
uation  constants  can  be  obtained  from  its  real  and  imaginary  parts. 

The  characteristic  impedance  is  calculated  as  follows.  After 
the  complex  propagation  constant  Is  obtained,  the  unknown  functions  f(x) 
and  g(x)  can  be  calculated  to  obtain  approximate  values  of  Ex  and  Jz> 
The  characteri stic  impedance  is  then  given  by: 
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I  E  (x)  dx 
a 


2  /  J  (x)  dx 


(  17  ) 

The  numerator  of  equation  (  17  )  represents  the  voltage  across  the  slot 
region  and  the  denominator  represents  the  total  current  flowing  on  the 
center  conductor  into  the  z  direction. 

COMPUTATIONAL  RESULTS 

According  to  the  above  theory,  a  FORTRAN  program  was  written  and 
the  integral  equation  (  14  )  was  successfully  solved.  A  CDC  Dual  Cyber 
170/370  at  The  University  of  Texas  Computation  Center  was  used  for  all 
calculations.  The  evaluation  of  the  singularity  of  the  matrix  Z  was  done 
by  looking  at  the  numerically  calculated  values  of  the  determinant  and 
the  condition  number  of  the  matrix.  These  values  were  obtained  by  the 
Gaussian  elimination  method  using  the  UNPACK  subroutines  in  The  Univer¬ 
sity  of  Texas  numerical  analysis  libraries.  The  accuracy  of  the  singu¬ 
larity  obtained  for  the  matrix  was  checked  using  the  singular-value 
decomposition  method  which  is  also  in  the  UNPACK  subroutines.  The 
search  for  the  complex  propagation  constant  was  done  by  using  the  complex 
root-search  routine  ZANLYT  in  the  IMSL  subroutines,  which  uses  the 
Muller's  method  with  deflation.  A  solution  was  usually  obtained  in  very 
small  number  of  iterations,  typically  6  to  10  times  depending  on  the 
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starting  point  of  the  iteration.  The  summations  in  the  integral  equation 
(  14  )  have  complicated  forms  and  therefore  need  to  be  truncated.  To 
check  the  accuracy,  the  behavior  of  the  results  with  respect  to  the  num¬ 
ber  of  terms  in  the  summations  over  n  and  the  number  of  matching  points 
was  investigated.  It  Is  not  easy  to  investigate  either  of  these  analyt¬ 
ically,  since  both  of  them  strongly  depend  on  the  dimensions  of  the 
structure.  Therefore,  the  convergence  of  the  solutions  was  checked 
numerically.  An  example  of  the  convergence  studies  is  shown  in  Fig.  3, 
which  is  calculated  for  a  typical  MIS  coplanar  waveguide.  The  figure 

shows  errors  relative  to  the  values  at  N  =30.  For  a  fixed  value  of  N  , 

P  P 

the  convergence  of  the  solutions  is  tested  with  respect  to  the  number  of 
terms  in  the  summations,  and  the  converged  values  are  plotted  in  the  fig¬ 
ure.  The  curves  are  not  very  smooth.  This  is  caused  by  the  nonuniform¬ 
ity  in  the  distribution  of  the  matching  points.  For  fixed  values  of 
material  constants,  such  as  dielectric  constants  and  conductivities,  the 
convergence  is  faster  if  the  thickness  of  the  insulated  layer  (  Region 
2  )  is  larger.  Thinner  insulated  layers  cause  the  fields  to  abruptly 
change  near  the  edges  of  the  conductors,  therefore  it  requires  more 
matching  points  and  more  terms  in  the  summations  to  fit  them.  If  there 
is  no  doped  semiconductor  layer  (  a  conventional  lossless  coplanar 
waveguide  ),  then  the  fields  change  very  smoothly.  Therefore,  It 
requires  rewer  matching  points  and  fewer  terms  in  the  summations  to 
obtain  precise  solutions. 

Computed  frequency  characteristics  of  the  slow-wave  factor, 
XQ/Xg,  (  XQ  Is  the  wavelength  in  free  space  and  X^  is  the  wavelength  in 
the  guide  )  and  the  attenuation  constant  of  the  MIS  coplanar  waveguide  is 


■Jib- 


error  relative  to  Np=30 


10 


o 


number  of  matching  points  Np 

Figure  3  Convergence  of  the  solutions 
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shown  In  Fig.  4.  Experimental  values  obtained  by  Hasegawa  are  also  shown 
for  comparison  [10].  They  are  in  reasonably  good  agreement.  As  the  fig¬ 
ure  shows,  slow-wave  propagation  exists  in  the  low  frequency  range,  and 
becomes  less  pronounced  as  the  frequency  increases.  The  transition  fre¬ 
quency  is  dependent  on  the  value  of  the  conductivity  of  the  doped  region. 
This  dependence  is  shown  in  Fig.  5.  There  exists  an  optimum  conductivity 
at  which  the  transition  frequency  of  the  slow-wave  region  extends  to  the 
largest  value  [12].  This  is  important  for  the  mm-wave  operation  of  the 
device. 

The  numerical  values  of  the  phase  and  attenuation  constants  are 
also  shown  in  Table  1  along  with  the  values  obtained  by  some  other  meth¬ 
ods  for  comparison  [15].  They  agree  well  one  another,  which  proves  the 
precision  of  the  proposed  method.  The  frequency  dependence  of  the  char¬ 
acteristic  impedance  is  shown  in  the  next  figure  (  Fig.  6  ).  The  values 
of  the  characteristic  impedances  agree  with  those  obtained  by  Wen  using  a 
quasi-TEM  result  for  the  lossless  case  [18].  The  different  tendencies  of 
the  two  curves  in  Fig.  6  will  be  discussed  in  Chapter  3  in  conjunction 
with  the  results  obtained  for  the  MIS  periodic  coplanar  waveguide. 


slow-wave  factor 


frequency  GHz 


Figure  4b  Attenuation  constant  versus  frequency 
o  :  experiment  [8] 
aB0.05mm,  b=0.5mm,  dj=0.4ym,  d2=3.0um, 
d 1 . Otwii •  e2*8*5eq«  e3—e^_13.0£q» 
a3=1.82xl04S/m 
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Table  1  Computed  phase  and  attenuation  constants 
obtained  by  several  different  methods 


Spectral -domain 

Mode-matchi ng 

Present  method 

f  =  0.1  GHz 

8/k0 

37.91 

35.67 

35.10 

a 

0.0466 

0.0386 

0.0322 

f  =  1.0  GHz 

e/ko 

31.61 

29.87 

28.81 

a 

1.76 

1.56 

1.51 

.1 

•1 

1 

'  1 
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imaginary  part  of  characteristic  impedance  Ohms 


CHAPTER  3  :  MIS  PERIODIC  COPLANAR  WAVEGUIDE 

As  shown  in  the  previous  chapter,  the  MIS  coplanar  waveguide  has 
a  very  large  slow-wave  factor  below  a  certain  critical  frequency.  For 
example,  the  result  shown  In  the  Fig.  4  indicates  that  the  electrical 
length  of  the  device  at  1  GHz  Is  about  30  times  larger  than  the  physical 
device  length.  However,  high  attenuation  may  make  the  application  of 
such  structures  difficult. 

In  this  chapter,  a  periodic  structure  Is  Introduced  to  reduce 
this  attenuation.  The  calculated  result  shows  that  the  MIS  periodic 
coplanar  waveguide  has  advantages  over  the  uniform  structure. 

APPROACH 

The  basic  theoretical  treatment  of  the  MIS  periodic  coplanar 
waveguide  consists  of  using  Floquet's  theorem  for  periodic  transmission 
lines.  To  do  this,  reasonably  accurate  values  of  both  propagation  con¬ 
stant  and  characteristic  Impedance  of  the  constituent  sections  of  each 
period  need  to  be  calculated.  Using  the  numerical  method  developed  in 
the  previous  chapter,  these  values  can  be  calculated  and  the  MIS  periodic 
coplanar  waveguide  can  be  analyzed  (  Fig.  1  ).  The  periodic  structure 
consists  of  two  different  sections:  one  is  an  MIS  coplanar  waveguide  and 
the  other  a  lossless  coplanar  waveguide  on  the  substrate  with  two  differ- 


ent  materials  (  Insulated  layer  and  semi-insulated  semiconductor 
substrate  ).  Therefore,  the  easiest  way  to  analyze  this  structure  Is  to 
treat  It  as  a  periodic  transmission  line  consisting  of  two  lines  with 
different  characteristics;  each  line  is  characterized  by  its  character¬ 
istic  impedance  and  electrical  length  (  Fig.  7  ).  In  this  approach,  the 
effect  of  the  geometrical  discontinuity  at  the  junction  of  two  sections 
is  neglected,  which  is,  however,  assumed  to  be  small. 

The  overall  propagation  characteristic  of  the  MIS  periodic 
coplanar  waveguide  is  approximately  calculated  by  applying  Floquet  s 
theorem  [19]: 


cose 


1  r  1  2 1 

=  cose,  cose,  —  —  +  —  sine. 

1  2  2Lz2  V 


sine. 


(  18  ) 


where 


6  =  2T1,  02  *2*2 
1  =  ll  *  i2  :  tV  lZ  1 ength  of  eacy  section 

The  characteristic  impedances  (  Zj,  Z2  )  and  the  propagation  constants 
(  JTj,  r2  )  of  two  sections  are  calculated  by  applying  the  numerical  meth¬ 
od  discussed.  The  propagation  constant  (  y  )  of  the  periodic  structure 
obtained  by  solving  eq.  (  18  )  is  again  a  complex  value,  from  which  we 
can  investigate  the  attenuation  constant  and  the  slow-wave  factor. 


Figure  7  Two  constituent  sections  of  the  MIS  periodic 

coplanar  waveguide  :  Each  section  is  characterized 
by  its  electrical  length  and  the  characteristic 
impedance. 


RESULTS 


The  computed  slow-wave  factors  and  attenuation  constants  for  MIS 
periodic  coplanar  waveguide  are  presented  in  Fig.  8  and  Fig.  9,  respec¬ 
tively.  Two  typical  cases  are  shown  in  the  figures.  A  remarkable  result 
is  obtained  for  the  case  03=10*  S/m  (  Fig.  8a  ),  where  the  extension  of 
the  frequency  range  of  the  slow-wave  propagation  is  observed.  Namely, 
the  slow-wave  factor  of  the  periodic  structure  becomes  greater  than  that 
of  uniform  MIS  coplanar  waveguide  at  frequencies  higher  than  10  GHz. 
This  crossover  of  the  slow-wave  factor  occurs  for  conductivities  larger 
than  a  certain  critical  value.  For  instance,  no  crossover  point  exists 
for  03=10*  S/m  (  Fig.  8b  ).  This  difference  can  be  explained  by  the 
behavior  of  the  doped  layer.  A  highly  conductive  doped  layer  tends  to 
act  like  a  conductor  at  high  frequencies  [2].  This  tendency  can  also  be 
seen  from  the  results  of  the  previous  chapter,  Fig.  6,  where  the  real 
part  of  the  characteri Stic  impedance  of  the  MIS  coplanar  waveguide 
becomes  smaller  at  high  frequencies  for  o^lO5  S/m.  Therefore  the 
effect  of  periodicity  becomes  strong  and  provides  slow-wave  propagation 
in  this  frequency  range.  In  fact,  if  the  doped  regions  of  the  MIS  peri¬ 
odic  coplanar  waveguide  are  replaced  by  conductor  strips,  the  resulting 
slow-wave  factor  follows  nearly  the  same  curve  at  high  frequencies.  On 
the  other  hand,  a  doped  layer  with  low  conductivity  acts  as  a  dielectric 
layer,  and  the  periodicity  becomes  weak  at  high  frequencies.  Therefore 
the  real  part  of  the  characteristic  impedance  of  the  MIS  coplanar  wave¬ 
guide  becomes  larger  at  high  frequencies  (  Fig.  6a  ).  In  this  case,  the 
periodic  structure  does  not  provide  the  enhancement  to  the  slow-wave 


frequency  GHz 


Figure  8a  Comparison  of  the  slow-wave  factor  of  the  MIS 
periodic  coplanar  waveguide  with  the  uniform 
MIS  coplanar  waveguide 
a=0.05mm,  b=0.5mm,  d^l.Outn,  d2=3,0ym, 

oo=105S/m,  Jl,  ={.o5:0. 1mm 


frequency  GHz 

Comparison  of  the  attenuation  constant  of  the 

MIS  periodic  coplanar  waveguide  with  the  uniform 

MIS  coplanar  waveguide 

a=0.05mm,  b= 0.5mm,  dj=1.0um,  d2=3.0pm, 

d^l.Omm,  E2=8.5e0,  €0=E^=13.Oeo, 

a,=105S/m,  2, =£o=0. 1mm 


behavior.  The  reduction  of  attenuation  constant  Is  also  more  pronounced 
for  the  case  03=10*  S/m  compared  to  the  other  case  (  Fig.  9  ).  There¬ 
fore,  the  MIS  periodic  coplanar  waveguide  is  useful  when  the  conductivi¬ 
ty  of  the  doped  regions  is  larger  than  a  certain  critical  value  around 
10*  S/m,  which  corresponds  to  the  doping  level  of  around  102*  /m3  in  gal¬ 
lium  arsenide  (  GaAs  )  substrate.  This  point  will  also  be  discussed  in 
Chapter  5. 


CHAPTER  4  :  EXPERIMENTS 


To  justify  our  theory,  a  simple  model  of  the  MIS  periodic  copla- 
nar  waveguide  was  fabricated  and  tested  in  the  frequency  range  of  40  MHz 
to  1  GHz.  Complex  propagation  constants  were  measured  using  a  conven¬ 
tional  Impedance  measurement  setup.  Instead  of  a  doped  semiconductor 
substrate,  graphite  powder  was  used  as  a  conductive  material.  The  gra¬ 
phite  powder  was  carefully  sandwiched  by  two  adhesive  plastic  sheets, 
and  placed  periodically  on  the  coplanar  waveguide  etched  on  a  circuit 
board,  which  has  approximately  the  same  dielectric  constant  as  the  plas¬ 
tic  sheets  (  Fig.  10  ) .  A  model  of  the  uniform  MIS  coplanar  waveguide  is 
also  fabricated  using  the  same  materials.  The  dimensions  and  other 
parameters  of  these  two  models  are  shown  in  Table  2. 

The  input  impedances  of  the  fabricated  lines  with  open  and  short 
terminations  were  measured  using  an  admittance  bridge  and  the  results 
were  converted  into  the  slow-wave  factors  and  attenuation  constants. 
The  conductivity  of  the  graphite  layer,  estimated  from  the  measured 
results  of  the  propagation  constants  of  the  model  of  the  uniform  MIS 
coplanar  waveguide,  was  about  20  S/m.  The  nominal  dielectric  constant 
of  the  graphite  layer  was  assumed  to  be  unity  since  the  graphite  is  a 
conductive  material.  The  experimental  results  are  shown  in  Fig.  11. 
Since  we  used  the  graphite  powder  to  create  conductive  layers  in  the 
coplanar  waveguide,  they  were  not  quite  uniform  in  thickness  and  densi¬ 
ty.  Nevertheless,  the  experimental  results  are  in  reasonably  good 
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Table  2  Dimensions  and  material  constants  of  the 
experimental  models 


common  parameters 

total  length  of  line 

212  mm 

width  of  center  conductor 

(  2a  ) 

2.0  mm 

width  of  slot 

(  b-a  ) 

4.0  mm 

thickness  of  substrate 

(  d3  ) 

3.0  mm 

dielectric  constants 

(  Ej,  e 2  and  e4  ) 

2.5  c 

0 

MIS  periodic  coplanar  waveguide 

thickness  of  insulator 

(  dj  ) 

0.13  mm 

thickness  of  graphite  layer  (  d£  ) 

0.02  mm 

length  of  each  section 

(  fcj  and  eg  ) 

8.0  mm 

number  of  periods 

12 

Uniform  MIS  coplanar  waveguide 
thickness  of  Insulator  (  d^  )  0.18  mm 

thickness  of  graphite  layer  (  d£  ) 


0.02  mm 


Figure  11a  Comparison  between  experiment  and  theory 
slow-wave  factor 
o  •  :  experimental  results 


MIS  coplanar  waveguide 


10  M 


100  M 
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Figure  lib  Comparison  between  experiment  and  theory 
attenuation  constant 


o  • 


:  experimental  results 


agreement  with  the  theoretical  curves.  The  crossover  of  the  slow-wave 
factor  does  not  occur  In  this  case  since  the  conductivity  of  the  graphite 
powder  is  too  small.  However,  this  experiment  verifies  the  theoretical 
calculations,  and  the  extension  of  the  slow-wave  propagation  range  is 
expected  if  appropriate  conductivity  is  obtained. 


CHAPTER  5  :  COPLANAR  SCHOTTKY  VARIABLE  PHASE  SHIFTER 


In  this  chapter,  discuss  a  Schottky-contact  coplanar  waveguide 
as  a  variable  phase  shifter  is  discussed.  This  device  can  be  constructed 
so  that  the  center  conductor  of  the  coplanar  waveguide  forms  a  Schottky 
contact  with  a  highly  doped  semiconductor  substrate.  Applying  a  proper 
dc  bias  voltage  forms  a  depletion  layer  at  the  contact,  which  has  a  much 
smaller  conductivity  than  the  surrounding  region.  This  situation  is 
similar  to  the  MIS  coplanar  waveguide,  and  slow-wave  propagation  is 
expected  on  such  a  line  (  Fig.  12  ).  Since  the  depletion  layer  thickness 
can  be  controlled  by  changing  the  dc  bias  voltage  on  the  center  conductor 
of  the  coplanar  waveguide,  there  is  the  possibility  of  making  a  electron¬ 
ically  variable  phase  shifter  using  this  configuration. 

The  theoretical  discussion  so  far  shows  that  the  MIS  coplanar 
waveguide  has  a  higher  attenuation  coefficient  at  high  frequencies.  The 
slow-wave  nature  also  disappears  as  the  operating  frequency  becomes 
higher.  However,  the  periodic  structure  proposed  in  the  previous  chap¬ 
ter  has  the  advantage  of  reducing  the  attnuation  constant  and  extending 
the  slow-wave  frequency  region.  The  idea  of  making  periodically  doped 
sections  on  the  semiconductor  substrate  can  also  be  used  in  the  Schottk¬ 
y-contact  coplanar  waveguide,  which  may  reduce  the  attenuation  and 
enhance  the  performance  of  the  variable  phase  shifter. 

In  the  present  chapter,  the  previously  developed  numerical  meth¬ 
od  will  be  used  to  determine  the  optimum  operating  conditions  for  both 


insulated  layer 


b  Schottky  coplanar  waveguide 
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uniform  and  periodic  Schottky-contact  coplanar  waveguides  as  electron¬ 
ically  variable  phase  shifters.  The  theory,  however,  can  not  be  directly 
applied  to  these  structures  since  they  have  curved  boundaries  at  the  bot¬ 
tom  of  the  depletion  layers.  In  order  to  be  able  to  utilize  the  theore¬ 
tical  method,  it  is  necessary  to  model  the  Schottky-contact  coplanar 
waveguide  by  the  MIS  coplanar  waveguide  with  the  insulated  layer  having 
the  same  hight  as  the  depletion  layer  of  the  Schottky  contact.  Since  the 
depletion  layer  is  very  thin  and  the  field  is  mostly  confined  under  the 
center  conductor,  the  error  introduced  by  this  model  Is  considered  to  be 
small.  Therefore,  the  performance  of  the  electronically  variable  phase 
shifter  can  be  studied  by  this  method. 


CONDITIONS  OF  OPTIMUM  OPERATION 


The  Important  parameters  to  be  considered  for  an  optimum  operat¬ 
ing  condition  are  the  dimensions  of  the  coplanar  waveguide,  the  conduc¬ 
tivity  and  the  thickness  of  the  doped  region,  and  the  depletion  layer 
thickness.  The  relation  between  the  dc  bias  voltage  V  and  the  thickness 
of  the  depletion  layer  dj  can  be  calculated  by  the  well  known  result  [5]: 


d 


1 


(  19  ) 


where  c  Is  a  permittivity,  q  the  unit  charge,  the  doping  concentration 
of  impurities,  V  applied  bias  voltage,  and  Vq  diffusion  potential  of  the 


Schottky  barier.  For  GaAs  substrate,  this  relation  is  plotted  in 
Fig.  13.  It  is  known  that,  at  a  fixed  frequency,  there  is  a  certain  val¬ 
ue  of  the  conductivity  of  the  doped  region  that  yields  a  maximum  range  of 
slow-wave  frequency  region  [12].  This  was  also  calculated  in  chapter  2. 
This  value  of  the  conductivity  also  gives  the  minimum  attenuation  over 
the  range  of  conductivities.  Therefore,  it  is  important  to  choose  an 
impurity  concentration  of  the  doped  region  so  that  the  conductivity  of 
this  region  becomes  this  optimum  value.  However,  the  impurity  level  is 
limited  to  a  certain  value,  because  the  breakdown  voltage  becomes  low 
when  the  impurity  concentration  is  large.  In  Fig.  13,  the  dotted  line 
indicates  the  critical  voltage  to  cause  the  breakdown  in  the  depletion 
layer,  which  is  obtained  by  using  a  known  value  of  breakdown  electric 
field  intensity,  Eg=4xl07  V/m  [20].  For  this  reason,  a  useful  value  of 
the  conductivity  of  the  doped  region  in  GaAs  substrate  is  limited  up  to 
10*  or  10s  S/m. 

Fig.  14  shows  the  behavior  of  the  optimum  conductivity  at 
10  GHz.  Three  curves  are  drawn  to  show  the  effect  of  the  thickness  of 
the  doped  region  d?.  It  is  interesting  that  Increased  value  of  d^  shifts 
the  point  of  minimum  attenuation  to  the  left,  i.e.,  the  optimum  conduc¬ 
tivity  is  lowered.  Also  the  minimum  value  of  the  attenuation  constant  is 
slightly  lowered.  Therefore  it  is  necessary  to  set  to  a  proper  value. 
This  tendency,  however,  disappears  when  d^  exceeds  a  certain  value, 
e.g.,  10  pm  in  this  case.  The  attenuation  constant  also  depends  on  the 
width  of  the  center  conductor  and  the  slot.  Fig.  15  shows  this  depend¬ 
ence.  The  relation  is  simple:  the  narrower  the  widths,  the  lower  the 
attenuation.  However  the  calculation  does  not  include  the  effect  of  con- 


depletion  layer  thickness 


Figure  13  Depletion  layer  thickness  versus  DC  bias 


current  for  GaAs  substrate  :  Dotted  line 
indicates  the  critical  voltage. 


attenuation  constant  dB/mm 
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conductivity  of  the  doped  region  S/m 


Figure  14  Behavior  of  the  optimum  conductivity 


f  =  10  GHz 


ductor  loss,  and  It  Is  expected  that  having  a  very  narrow  center  conduc¬ 
tor  causes  large  attenuation.  Also,  there  is  a  lower  limit  on  the  slot 
width  because  the  breakdown  voltage  of  the  semiconductor  substrate 
becomes  smaller  for  narrower  slots. 

There  is  one  more  parameter  to  consider:  depletion  layer  thick¬ 
ness  dj.  The  thickness  d^  is  determined  by  the  dc  bias  voltage,  and  as 
was  discussed  earlier,  it  changes  the  phase  shift  of  the  device.  Howev¬ 
er,  it  also  changes  the  attenuation  constant.  This  is  inherent  to  this 
structure,  but  it  is  desirable  to  keep  this  change  as  small  as  possible. 
Fig.  16  is  drawn  for  2  different  values  of  dj.  The  figure  shows  that  the 
small  values  of  d^  yields  large  attenuation.  Also,  the  point  of  minimum 
attenuation  moves  as  d^  changes.  This  means  that  we  would  like  to  use 
point  A  in  the  figure  rather  than  B  to  keep  the  total  change  in  the 
attenuation  small . 

RESULTS 

According  to  the  above  discussion,  a  model  can  be  constructed 
for  the  optimum  operation  of  a  Schottky  variable  phase  shifter.  A  GaAs 
substrate  is  chosen  because  of  its  large  electron  mobility  and  breakdown 
electric  field  intensity.  Various  constants  for  GaAs  substrate  used  in 
the  calculation  are  shown  In  Table  3.  The  width  of  the  center  conductor 
of  the  coplanar  waveguide  is  chosen  to  be  2  pm  (  a=l  um  ).  In  this  case, 
the  doping  density  in  the  doped  region  Is  N^SxlO23  /m3,  which  yields  a 
conductivity  of  02=4.1x10*  S/m.  This  gives  the  lowest  attenuation  and  a 


attenuation  constant  dB/mm 
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conductivity  of  the  doped  region  S/m 


Figure  16  Attenuation  constant  versus  conductivity 
of  the  doped  region  :  Point  A  is  better 
than  B  since  the  attenuation  varies  less, 
f  =  10  GHz 
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Table  3  Various  constants  for  GaAs  substrate 
used  In  the  calculation 

dielectric  constant  13.0 

mobility  8500  cm*/V-sec 

breakdown  field  intensity  4.0  *  107  V/m 

impurity  density  (doped  region)  3.0  x  10**  /m* 


depth  of  doped  region  (dg) 


10  ym 


calculated  maximum  depletion  layer  thickness  of  djmax=0.12  ym  when  the 
dc  bias  voltage  is  V+Vq=1.9  V.  This  Is  the  maximum  voltage  that  can  be 
applied  In  this  situation  (  breakdown  voltage  ). 

Using  these  values,  the  required  length  of  the  device  and  the 
total  attenuation  in  the  device  are  calculated.  The  depth  of  the 
depletion  layer  is  varied  from  0.03  to  0.09  ym,  which  corresponds  to 
adjusting  V+Vq  from  0.2  to  1.7  V.  The  device  length  is  determined  so 
that  the  difference  in  phase  of  the  output  signal  becomes  180  degrees 
when  the  dc  bias  voltage  Is  changed  from  minimum  to  maximum  value  (  0.2 
to  1.7  V  ).  The  results  are  shown  in  Fig.  17.  This  figure  Is  drawn  as 
follows;  for  each  value  of  b,  the  total  device  length  required  for  180 
degree  phase  shift  is  calculated  based  on  the  propagation  constant 
obtained  by  our  theoretical  calculation.  The  total  attenuation  In  the 


device  Is  then  obtained  by  multiplying  the  maximum  attenuation  constant 
for  the  operating  condition  by  the  device  length.  The  curves  show  that 
when  the  slot  width  becomes  narrower,  the  total  attenuation  becomes 
smaller  in  spite  of  the  longer  device  length  required.  The  figure  also 
shows  that  the  total  attenuation  is  smaller  for  the  periodic  structure 
compared  to  the  uniform  device.  However,  this  advantage  becomes  less 
when  the  slot  width  becomes  very  narrow.  The  performance  of  the  periodic 
structure  could  have  been  much  better  if  the  conductivity  of  the  doped 
region  were  higher.  We  could  not  do  this  because  of  the  low  breakdown 
voltage  of  GaAs.  Also,  because  of  the  same  limitation,  we  can  not  meet 
the  last  condition  given  in  the  previous  section:  minimize  the  variation 
of  the  attenuation  constant  over  the  operating  dc  bias  voltage.  There¬ 
fore,  if  a  better  substrate  material  were  used  (  i.e.  one  with  high  elec- 


Figure  17a  Required  device  length  for  various  slot  width 
a  *  1.0  yin,  Parameters  for  GaAs  substrate  are 
shown  in  Table  1-3. 


tron  mobility  and  breakdown  electric  field  ),  an  Improvement  of  the 
performance  of  the  periodic  variable  phase  shifter  would  be  possible. 

Finally,  the  total  attenuation  of  the  device  Is  plotted  as  a 
function  of  the  frequency  In  Fig.  18.  For  a  fixed  frequency,  the 
required  device  length  Is  calculated  and  the  total  attenuation  is  then 
computed  for  that  length.  At  40  GHz,  the  calculated  total  attenuation  is 
1.27  dB  for  uniform  device  and  1.20  dB  for  the  periodic  device.  Although 
these  figures  do  not  Include  conductor  loss,  they  indicate  quite  satis¬ 
factory  performance  of  these  devices  at  millimeter-wave  frequencies. 


Figure  18  Total  attenuation  versus  frequency 
a  =  1.0  um,  b  =  1.5  ym 
Parameters  for  GaAs  substrate  are  shown 


In  Table  1-3. 


CHAPTER  6  :  CONCLUSIONS 


In  Part  I,  the  theoretical  development  of  the  coplanar,  $•  hottky 
variable  phase  shifter  has  been  accomplished.  An  efficient  computer 
program  was  developed  to  calculate  the  trasmissi on-line  parameters  of 
the  MIS  slow-wave  coplanar  waveguide.  In  order  to  reduce  the  inherent 
loss  in  such  a  device,  a  periodic  structure  was  proposed  and  analyzed  by 
conventional  transmission-line  theory.  A  simple  experiment  verified  the 
applicability  of  this  method.  Extensive  calculations  were  performed  to 
find  the  optimum  conditions  of  the  phase  shifters  based  on  the  GaAs  subs¬ 
trate.  Calculations  show  that  optimized  180-degree  phase  shifter  can  be 
operated  at  40  GHz  with  a  loss  as  small  as  1.2  dB. 


APPENDIX  A  :  DERIVATION  OF  THE  INTEGRAL  EQUATION 


In  Appendix  A,  a  detailed  derivation  of  the  Integral  equation 
(  14  )  will  be  given.  The  potentials  41  and  0  In  each  region  are  given  by 
eq.  (  7  )  through  (  11  ).  Using  these  expressions,  we  can  satisfy  all 
the  boundary  conditions  to  obtain  the  Integral  equation  (  14  ).  To  do 
so,  we  need  to  utilize  eqs.  (  5  )  and  (  6  ),  which  can  be  rewritten  as 
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The  boundary  conditions  to  be  satisfied  is  the  continuity  of  the  x  and  z 
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From  these  relations,  we  obtain 
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where 
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Using  (A-3),  (A-5),  and  (A-7),  we  obtain 
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Therefore,  we  can  obtain  the  ratio  of  C„  and  C  1  as  follows: 

n  n 
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These  are  the  coefficients  appearing  in  eq.  (  12  ). 

Next,  by  applying  the  boundary  conditions  on  Ex  and  Ez  at  y=0, 
the  following  relations  can  be  easily  obtained. 
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Therefore  the  potentials  In  Region  2  can  be  written  as  in  eq.  (  12  ). 


In  order  to  obtain  the  Integral  equation  (  14  ),  we  use  the 
remaining  boundary  conditions  In  (  13  )  and  (  15  ).  Using  the  condition 
on  Exj,  we  obtain 
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which,  by  the  Fourier  analysis,  becomes 


1  2  fb 

-  An  “in  pn  ‘  jjr  B"  =  “  >  f(x)  s1n(J"X  dx 


Jut 


1 


w  a 


(A— 9) 


Similarly,  using  the  condition  on  Hx, 
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Using  these  results,  we  can  obtain  the  coefficients  An  and  Bn: 
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We  now  substitute  An  and  Bn  Into  the  rest  of  the  boundary  conditions: 
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These  relations  can  be  rewritten  as  the  Integral  equation  (  14  )  using 
the  Fourier  analysis: 

Z  C  pn  /  f(x')  s1n3nx'  dx'  +  Qn  /  g(x')  cos^x'  dx'  ]  cos3nx  =  0 
ns  c 

0<x<a,  b<x<w 

Z  [  Rn  /  f(x')  s1n3nx'  dx'  +  Pp  /  g(x')  cos3nx'  dx'  ]  s1n&nx  =  0 
ns  c 

a<x<b 

(1-14) 

where 


APPENDIX  B  :  PROGRAM  LISTING 


Appendix  B  contains  the  FORTRAN  program  code  for  the  analysis  of 
the  propagation  characteristics  of  the  uniform  and  periodic  MIS  coplanar 
waveguide.  The  subprograms  CGECO,  CGEDI,  CGESL,  and  CSVDC  are  found  in 
the  UNPACK  library,  and  ZANLYT  is  found  in  the  IMSL  library  in  The  Uni¬ 
versity  of  Texas  Numerical  Analysis  Library. 
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PROGRAM  CPWINT(T?Y, OUTPUT, TAPEl. TAPED=TTY, TAPEo=TTY) 


COMPUTATION  CiF  THE  COMPLEX  PROPAGATION  CONSTANT 
OP  THE  UNIFORM  AND  PERIODIC  MIS 
COPLANAR  WAVEGUIDE 


STRUCTURE  : 

THE  RIGHT  HALF  OF  THE  CROSS  SECTION  OF  THE 
STRUCTURE  IS  SHOWN  BELOW 


M.  WALL  <  CENTER  ? 

EPS1  (AIR) 


D1 


D2 

D3 


<  - W - > 

<  - B - > 

<  -  A  -  > 


EP52  ( INSULATOR  > 


EPS3.SIG3  (DOPED  REGION) 


EPS4  (SEMI -INSULATED) 


EPSS  (AIR) 


WALL 


A  ...  HALF  WIDTH  OF  CENTER  STRIP 
B  ...  A  +  WIDTH  OF  SLOT 
W  ...  HALF  WIDTH  OF  OUTER  GUIDE 
D1  .  .  .  THICKNESS  OF  INSULATOR 
D2  .  .  .  THICKNESS  OF  DOPED  REGION 
D3  THICKNESS  OF  SI  SUBSTRATE 

MM AX  . . .  NUM3ER  OF  TERMS  IN  SUMMATIONS 
NP  .  .  NUMBER  OF  POINTS  TO  EVALUATE 

INTEGRAL 

SOME  VARIABLES  ARE  PRESET  BY  THE  PROGRAM  TO 
CERTAIN  VALUES  ( GAAS  IS  ASSUMED).  WHEN  YOU 
START  CALCULATION  FOR  A  DIFFERENT  STRUCTURE, 
CONVERGENCE  OF  THE  SOLUTIONS  SHOULD  BE  INVESTI¬ 
GATED  BY  CHANGING  THE  VALUES  OF  MM AX  AND  NP. 

(  MM  AX  .  LE.  1200,  NP  LT.  50  ) 

THE  SOLUTION  IS  OBTAINED  BY  EVALUATING  THE 
DETERMINANT  OF  THE  MATRI *  Z  IN  FUNCTION  SUB¬ 
PROGRAM  DETERM.  THE  ACCURACY  OF  THIS  METHOD 
CAN  BE  CHECKED  BY  CHOOSING  THE  OPTION  7  AFTER 
THE  CALCULATION  OF  THE  PROPAGATION  CONSTANT 
IT  USES  THE  SINGULAR  VALUE  DECOMPOSITION 
TO  ACCURATELY'  EVALUATE  THE  CONDITION  NUMBER  OF 
THE  MATRIX  3  CONDITION  NUMBERS  ARE  CALCULATED 
Bv  THE  PROGRAM  (RCOND) .  THE  1ST  ONE  SHOULD  BE 


o  n  r»  o  ono  ooooooooo 
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MUCH  SMALLER  THAN  THE  REST  OF  THEM 


LIBRARIES  : 

AALIB  ( CGECO.  CGEDI. CCESL,  CTAN, PHASE, 
CSVDC. CSROT) 

IMSLIBF  (IANLYT) 


COMMON  /BLOCKO/A,  B-  W,  Dl.  02,  D3,  AKO,  PI,  ZMU,  ZEP,  OMEGA 
COMMON  /BLOCKl/EPSl, EPS2. ZEPS3, EPS4, EPS5 

COMMON  /BL0CV\2 /SETA  (1200) 

COMMON  /BLQCK3/MMAX,  NP.  Jl,  J2,  J3,  I  IMP,  I  FACT,  ZIMP 
COMMON  /BLQCK4 / X ( SO ) , XP  <  5C ) 

COMPLEX  ZEPS3.  Z J, DETERM 
COMPLEX  ZDET , ZC (2), Z IMP. C < i 2) , PCONST ( 2i 
REAL  L<2>,T0<4) 

DIMENSION  INFER •2) 

EXTERNAL  DETERM 

BASIC  CONSTANTS. 

ZJ  -  <0.  0.  1.  0) 

PI  -  3.  1415926S35897 
ZMU  =  4.  OE-7  *  PI 
ZEP  »  8.  854E-12 
DBN  *  8.  6S59E-3 

I FLAG 1  >  0 
I FLAGS  «  0 
IFLAG3  -  0 
I SEC  TN  -  1 
EPSl  -  1.0 
EPS2  -  13.  O 
EPS3  -  13.  O 
EPS4  -  13.0 
EP35  =1.0 
MMAX  =  500 
NP  =  25 
0  *  2.  0 
DO  1  1=1,4 
TG< I )  =  O  0 
1  CONTINUE 

ZC(1 )  ■  <0.  0,  0.  0) 

ZC<2)  ■  (0.  0,  0.  C) 

SIG3  =  0.  0 
FREQ  =  0.  0 
MM IN  =  0 
GAM  =00 
ATTEN  ■  0.  0 
SHIFT  ■  0.  0 

MAIN  CONTROL 

1000  WRITE  <6. 100) 

100  FORMAT  (  5/.5X,  '-*■****  MAIN  MENU  ***•*' 

1  //5X. '1  ...  CALCULATION  OF  MIS  COPLANAfi  WAVEGUIDE' 

2  /SX;  '2  .  CALCULATION  0"  PERIODIC  STRUCTURE' 

3  /SX,  '10  S'J'OP'  ) 

READ  <5,  ■*  )  1CTPL 

IF  ,  1CTRL  EO  1  >  50  TO  ZCOO 

IF  (  ICTRL  EQ  2  )  50  TO  -if  OO 
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IF  (  1CTRL  EG  10  l  STOP 
WRITE  (t>,  110; 

110  FORMAT  <  ) 

GO  TO  1000 
C 

C  CALCULATION  OF  MIS  COPLANAR  WAVEGUIDE. 

C 

2000  IF  1  I FLAG 1  EQ  0  >  GO  TO  2100 

WRITE  (6.120)  MnAX.  NP,  G.  Ul,  J2,  J3,  MMIN. 

1  EPS).£P=2.  EP53,  EPS4,  EPS5, 

2  A.  B.  W,  Dl.  D2,  D3,  SIG3.  FREQ 

120  FORMAT  <  5/.5X.  '»**♦**  MIS  COPLANAR  WAVEGUIDE  ******' 

1  /  /  5  X .  'MMAX  =  14.  '.  NP  *  '.  2,  Q  »  ',  F4.  1 

2  /5X  'DISTRIBUTION  '.3(X. I2>.'  (MMIN  14, ' ) ' 

3  /SX, 'DIELECTIC  CONSTANTS  :  '.F4  1,4(2X,F4. 1) 

4  /SX.  'A  *',E9  3.  '.  S  « ER  3.  ' .  W  »',E9.  3,  '  <M>  ' 

5  /SX.  'Dl  «'.E9  3.  D2  E9.  3,  D3  =',  E9.  3,  '  (M)  ' 

6  /5X,  'SIG3  3,  '.  FREG  «  F6.  2.  '  (GHZ)'  ) 

C 

IF  <  IFLAG2  EG  0  )  GO  TO  2010 
WRITE  (6.122)  GAM.  ATTEN 

122  FQRhAT  (  ,'5X,  'OAM/KO  -  ',F9  5.  '  ATTEN  =  '.E11.5, 

1  ’  DB/MM'  ) 

IF  (  IFLA02  EG  1  )  GO  TO  2010 
WRITE  (6.124)  ZIMP 

124  FORMAT  (  /SX, 'CHARACTERISTIC  IMPEDANCE  •  '.E11.5, 

1  2X,  Ell.  S  ) 

2010  WRITE  (6.  12b, 

126  FORMAT  (  /SX.  '1  ...  CALCULATE  PROPAGATION  CDNST.  ' 

1  /SX, '2  ..  CALCULATE  PROP.  CONST.  AND  IMPEDANCE' 

2  /SX. '3  .  .  CALCULATE  CHARACTERISTIC  IMPEDANCE' 

3  /SX,  '4  ...  CHANGE  1  PARAMETER ' 

4  /SX. '5  ...  CHANGE  DIMENSIONS,  SI03.  AND  FREG' 

5  /SX, '6  ...  CHANGE  all  PARAMETERS' 

6  /SX.  '7  ...  CHECK  WITH  CSvDC ' 

7  /SX,  '10  .  .  CO  BACK  TO  MAIN  MENU'  ) 

2020  READ  (S.  *  )  ICTRL 

IF  (  ICTRL  EG  1  .  OR.  ICTRL.  EG.  2  )  GO  TO  2500 

IF  (  ICTRL  EG  3  >  GO  TO  2030 

IF  (  ICTRL  EG.  4  )  CO  TO  2400 

IF  (  ICTRL  EG  5  )  GO  TO  2100 

IF  (  ICTRL  EG.  6  )  GO  TO  2200 

IF  (  ICTRL  EG.  7  )  GO  TO  3500 

IF  (  ICTRL  .EG  10  )  GD  TO  1000 

WRITE  (6, 110) 

GO  TO  2020 
C 

2030  IF  (  I FLAGS  .EG  1  )  GO  TO  3300 
WRITE  (6. 110) 

GO  TO  2020 
C 

2100  WRITE  (6,  130) 

130  FORMAT  (  5/, SX.  ******  INPUT  PARAMETERS  ******'  ) 

CO  TO  2300 
2200  WRITE  (6. 130) 

WRITE  (6.  140) 

140  FORMAT  (  /SX, 'MMAX,  NP,  Q'  ) 

READ  (5.  *  >  MMAX.NP.G 
WRITE  (6. ISO) 

150  FORMAT  (  /SX.  DIELECTRIC  CONSTANTS  1-3'  ■ 

READ  1 5,  *  >  EF-21.  EP22.EPS3.  EPS4.  EPS5 

2300  WRITE  <6,  160) 
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FORMAT  <  /5X,  'A  B.  W  <M>  '  : 

pead  <5,  *  >  a.  e.w 

WRITE  (6.  1TO) 

FORMAT  <  /5X. '01.  02.  03  <M' '  ) 

READ  <5,  *  )  01 .  02.  03 

WRITE  <i, ISO) 

FORMAT  <  /5X, 'SI53  (S'M),  FREQ  (CHZ)'  ) 

READ  (5,  *  )  S!G3. FREQ 

I FLAG 1  -  1 
IFLAG2  «  0 

CALL  POINTS  <X,UP.  Jl.  02.  J3.  Q) 

MM IN  -  W/(X< Ul+2'-X( Jl+1 ) >+l 
00  TO  2000 

IFLAG2  -  0 
WRITE  (6.  190) 

FORMAT  <  /5X, 'INPUT  NAME  OF  THE  PARAMETER.  '  ) 
READ  (5.400)  NAME 
FORMAT  (  A4  ) 

WRITE  (6.200) 

FORMAT  (  /5X. 'ITS  VALUE'  ) 

READ  (5,  *  )  VALUE 


IF 

( 

NAME 

EG 

4HMMAX 

) 

MMAX  »  VALUE 

IF 

( 

NAME 

EG. 

4HNP 

) 

GO  TO  2410 

IF 

( 

NAME 

EQ. 

4HG 

) 

G  -  VALUE 

IF 

( 

NAME 

.  EQ. 

AHA 

) 

A  =  VALUE 

IF 

( 

NAME 

.  EQ. 

4HB 

) 

B  <*  VALUE 

IF 

( 

NAME 

EQ. 

4HW 

) 

W  «  VALUE 

IF 

( 

NAME 

.  EQ. 

4HD1 

) 

D1  -  VALUE 

IF 

( 

NAME 

EQ. 

4HB2 

) 

D2  -  VALUE 

IF 

( 

NAME 

.  EQ. 

4HD3 

) 

D3  »  VALUE 

IF 

( 

NAME 

EQ. 

4HEPS1 

) 

EPS!  •<  VALUE 

IF 

( 

NAME 

.  EQ. 

4HEF32 

) 

EPS2  -  VALUE 

IF 

( 

NAME 

.  EQ. 

4HEPS3 

) 

EP53  -  VALUE 

IF 

< 

NAME 

EQ. 

4HEPS4 

) 

EPS4  -  VALUE 

IF 

< 

NAME 

EG. 

4HEPS5 

) 

EPS 5  -  VALUE 

IF 

( 

NAME 

.  EG. 

4H3IG3 

) 

SI 03  -  VALUE 

IF 

( 

NAME 

.  EG. 

4HFREG 

) 

FREQ  -  VALUE 

.  J3,  Q) 
X)  )  +  l 


GO  TO  2000 
2410  NP  *  VALUE 

CALL  POINTS  <X,NP. J1.J2 
MM IN  -  W/<X(Ul+2)-X(Ul+ 
GO  TO  2000 

r* 

:  START  ROOT  SEARCH. 


IF  (  ICTRL  EQ  i  )  IFLAG2  ■  1 
IF  <  ICTRL  EG  2  )  IFLAG2  -  2 
A*0  «  20.  0*PI*FRE(i/2.  998 
OMEGA  »  2.  0E9»Pl*FREG 

ZEP53  «  EPS3  -  (0  0.  J  0  >*SIG3/ < OMEGA* ZEP ) 
DO  5  1-1 ,  NP 

XP(I)  -  (X(l:+Xi t  +  l)>/2. 0 
CONTINUE 

DO  10  N-l.MMAX 

BETA(N)  »  <FL0A*'N>-0  5)*PI/W 
CONTINUE 


uuu 
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M3 I G  «  4 
UMAX  =  20 
WRITE  <6,22 0) 

220  FORMAT  <  /5X, 'INPUT  INITIAL  GUESS'  GAM/kO, 

1  ATTEN  (DB/MM)  > 

READ  <5,  *  )  GA:«,  ATTEN 

CXU)  -  GAM*aKO-ZJ*ATTEN/DBN 
2800  WRITE  (6,230) 

230  FORMAT  <  / 5X,  ROOT  SEARCH  STARTED. . .  ' 

1  //5X.  GAM/KO-'.  6X,  '-ATTEN-  ,  11X,  '-DET-', 

2  14X,  '-RCOND- '  ) 

C 

CALL  ZANLYT  (DETERM,  0  0,  NSIG.  0.  1 ,  1 .  CX.  I TMAX,  INFER.  IER  > 
C 

IF  (  REALiCX(l)  .  LT.  0.0  )  CX(1>  -  -CX(1) 

GAM  =  REAL ( CX <  * : ■  /  AKO 
ATTEN  =  -  A I MAG • C  X ( 1 > ) »DBN 
C 

IF  <  IER  .  EG  0  /  GO  TO  3100 
WRITE  (6,250)  GAM,  ATTEN,  IER 

250  FORMAT  (  //5X.  '•►*****  TERMINAL  ERROR  (ZANLYT)  ****•*•*' 

1  /5X,  'GAM/KO  =  F9.  5,  ',  ATTEN  «  ',£11.  5, 

2  /5X. '  DB/MM,  IER  «  ■ . 13 

3  ,'SX,  '1  ...  CONTINUE  CALCULATION' 

4  ,'5X,  '2  ...  D:-FERENT  INITIAL  GUESS' 

5  /5X, '3  ...  CHANGE  PARAMETERS'  ) 

2900  READ  <5,  *  )  IC^RL 

IF  (  ICTRL  .  EG  )  GO  TO  3000 

IF  <  ICTRL  . EG  2  )  GO  TO  2700 

IF  (  ICTRL  EG  3  )  GO  TO  5000 

WRITE  (6.110) 

GO  TO  2900 
C 

3000  WRITE  (6.260) 

260  FORMAT  (  /5X, ' I7MAX '  ) 

READ  <5.  *  )  I7MAX 
GO  TO  2800 
C 

3100  WRITE  (1,500)  MMAX. NR, G. 

1  EPS1 ,  EPS2,  EPS3,  EPS4,  EPS5, 

2  A,  B.  W,  Dl,  D2,  D3,  SIG3,  FREQ. 

3  GAM, ATTEN 

500  FORMAT  (  4/,  10X.  'MMAX  *  ',14,  NP  «■  ',12,  ',  G  m  ', 

1  F4  1/10X,  'DIELECTRIC  CONSTANTS  ' , F4.  1 , 4 < 4X. F4.  1 > 

2  /  10X ,  'A  -  ' ,  E9.  3 ,  ' ,  B  *'.  E9.  3,  ',  W  »',E9.3 

3  /10X,  'Dl  «',E9  3.  D2-',E9.  3,  ',  D3-',E9.  3 

4  /10X,  'SIG3  «',E9  3.  ',  FRED  -  ',F6  2.  '  GHZ' 

5  //10X,  'GAM/KO  *  ',  FlO.  6,  ',  ATTEN  -  '.£12.  6, 

6  '  DB/MM'  ) 

C 

IF  (  IFLAG3  . NE  0  )  GO  TO  3300 

IF  <  IFLAG2  EG  1  )  GO  TO  2000 

CALCULATE  IMPEDANCE 

3300  IFLAG2  -  2 
I  IMP  »  1 

C  X i 2 )  *  DETERM  C  X ( 1 )  , 

PCONST  ( ISECTN  >  =  CX(H 
TG (2-* (I SEC TN-1  )•  l  )  =  GAM 
TGO<  ISECTn-l  >-.:•)  =  ATTEN 
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ZC<  ISECTN)  =  ZJfM 
WRITE  <1.510  CIMP 

510  FORMAT  <  /10X.  'CHARACTERISTIC  IMPEDANCE  - 
1  £12.  o,  2X-  El?  6  ■ 

C 

IF  l  IFLAG3  EQ  1  )  GO  TO  4000 
00  TO  2000 
C 

C  CHECK  ACCURACY  BY  SINGULAR  VALUE  DECOMPOSITION. 

C 

3500  IF  <  IFLA02  NE  0  )  00  TO  3510 
WRITE  (6,110) 

00  TO  2020 
C 

3510  I  IMP  -  10 

CX<2>  =  DETERM  <CX(1> > 

C X  f  2 >  -  DETERM  • C X < 1 > * < 1.  O+lO  0** < -NSIG > ) > 

C  X  <  2 )  «  DETERM  ■-  C  X  <  1  •  *  1 1.  O-lO.  0**  (-NSIC ;  )  i 
CO  TO  2010 
C 

C  MIS  PERIODIC  COPLANAR  WAVEGUIDE, 

C 

4000  WRITE  (6,290)  7G < 1 > . TG<2> . ZC < 1 > , TG<3> . TC < 4 ) , ZC <2> 
290  FORMAT  <  5/,  5X.  -»#■*•#**  PERIODIC  STRUCTURE  ■******' 

1  //5X,  'SECTION'.  5X,  'GAM'.SX.  'ATTEN',  14X.  ZC' 

2  /ex,  '1  SX.  F9.  5.  2X.  Ell.  5.  2X,  E9.  3,  X.  E9.  3 

3  /8X,  '2',  5X,  F9.  5,  2X.  Ell.  5.  2x,  E9  3,  X.  E9.  3  ) 

IF  (  IFLAG3  NE  2  )  GO  TO  4010 

WRITE  (6,292)  GAM, SHIFT,  ATTEN 
292  FORMAT  (  /5X,  'CAM/KO  »  '.F9  3,  '  ('-F10.  5,  '  > ' 

1  /5*,  'ATTEN  *  ,  Ell.  5  ) 

4010  WRITE  (6. 294) 

294  FORMAT  <  /SX. '1  CALCULATION  FOR  SECTION  1' 

1  /SX,  '2  ...  CALCULATION  FOR  SECTION  2' 

2  /5X,  '3  ...  CALCULATE  PROPAGATION  CONST.  ' 

3  /5X, '10  .  GO  BACK  TO  MAIN  MENU'  ) 


4100  READ 

(5,  * 

)  IC 

TRL 

IFLAG3  -  1 

IF 

< 

1CTRL 

.  EQ 

l  ) 

GO 

TO  4200 

IF 

< 

ICTRL 

.  EQ 

2  ) 

GO 

TO  4300 

IF 

< 

ICTRL 

EQ 

3  ) 

GO 

TO  4500 

IF 

( 

ICTRL 

.  EQ 

10  ) 

GC 

1  TO  4400 

WRITE  (6,110) 

GO  TO  4100 
C 

4200  ISECTN  *  1 
GO  TO  2000 
4300  ISECTN  -  2 
GO  TO  2000 
4400  ISECTN  -  1 
IFLAG3  «  0 
GO  TO  1000 
C 

4300  IFLAG3  -  2 

WRITE  (6,300) 

300  FORMAT  (  / /5X ■ ‘LENGTH  OF  EACH  SECTION  (Mi'  ) 
READ  (5,  *  )  L‘ : ) , L (2) 

C 

CALL  PERIOD  (AKC  L, ZC, PC0NS7. GAM,  SHIFT, ATTEN) 

C 

WRITE  (1,520)  l.<  1 ),  L<2),  GAM.  SHIFT- ATTEN 

520  FORMAT  (  /10X,  'PERIODIC  STRUCTURE  (  LENGTHS  .  ' 


uouu o  uouuu 
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1  E12.  6.  4X.EJ2  o/ '10X,  'GAM/KO  - 

2  r  10.  5.  '),  ATTEN  -  ',£12.  i>  > 
GO  TO  4000 

END 


COMPLEX  FUNCTION  DETERM  (G> 

CONSTRUCTS  MATRIX  AND  CALCULATES  ITS  DETERMINANT. 
ALSO  SOLVES  FOR  ELECTRIC  CIELD  AND  CURRENT 
DENSITY  TO  CALCULATE  CHARACTERISTIC  IMPEDANCE. 

COMMON  /3L0CK0/A.  E.  W,  Dl,  D2.  D3,  AKO,  PI,  IMU,  ZEP,  OMECA 
COMMON  /BLQCK1/EPS1, EPS2,  ZEPS3, EPS4, EPBS 
COMMON  /SLOCkE/BETA! 1200/ 

COMMON  /BL0CK3/MMAX,  NP,  J) ,  J2,  J3. IIMP. IFACT , Z IMP 
COMMON  /BL0CK4/X  !  SO ) , XP  <  5C  > 

COMPLEX  G, GAMMA. GAMMA2, ZEPS3, ZIMP, CTAN 
COMPLEX  ZOO.  30).  WORK  <  50) .  SC  (  50 ) ,  DET  ( 2 ) 

COMPLEX  ALPHAO),  0.  P,  0.  R.  3.  T.  U,  V.  PCN,  PDN,  SG3ETA 
COMPLEX  PI,  01,  fil,  R.2.  R3,  R4.  R5.  R6,  SI,  S2 
INTEGER  IPUT  <  50  > 

REAL.  EA!  50/ .  EP  <  50) 

LDZ  50 

DBN  -  8.  68S9E-3 

GAMMA  »  G 

GAMMA2  *  GAMMA*GAMMA 
AK02  =  AKO*AKO 

DO  10  1  =  1,  NP 
DO  10  J=l, NP 
10  Z(U.  I )  «  (0.  0.  0  0) 

DO  200  N=1,MMAX 
SQBETA  =  BETA! N‘ *BETA<N) 


P  =  GAMMA2  -r  SOiETA 
ALPHA  < 1 )  =  CSQRTiP-EPSl *AK02 1 
ALPHA ( 2 )  »  CSQRT ;EFS2*AK02-r ; 
ALPHAO)  =  CSQRT  i  ZEPS3*AK02-s:  > 

ALPHA ! 4 )  =  CSQRT ; EPS4*AK02-P  > 

ALPHA! 5)  =  CSQRT (P-EPS5*AK02> 

0  =  CTAN (ALPHA! 4, *D3> 

PI  =  ALPHA ! 4 ) /ER S 4  *  0*ALPHA ! 5 ) /EPS5 
01  =  0* ALPHA! 4 ) / EPS4  -  ALPHA < 5 ) /EPS5 
R1  =  EPS2/ZEPS3-*.4LFHA  C  3 )  /  ALPHA  ( 2 ! 

R2  ZEPS3/EP34-*  ALPHA < 4 )  / ALFHA ( 3 ) 

R3  =  CTAN< ALPHA! 2 >*D1 ) 

R4  =  CTAN! ALPHA  3)*D2) 

R5  *  — R1*R3 
R6  =  -R2*R4 

51  •=  1 .  0+Q1/P1*!3* 

52  =  R4+Q1 /P 1*R2 

PCN  =  -!S1+S2*R3  /<S1*R3+S2*R1 > 

PI  «=  ALPHA! 4)  +  0*ALPWA(5> 

01  =  0* ALPHA' 4)  -  ALPHA! 5) 

PI  =  ALPHA  < 3 i / A_PHA ! 2  > 

R2  -  ALPHA ! 4 ! /ALPHA ( 3 ) 
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R5  *  -R1#R3 

R6  »  -R2*R4 

51  =  1.  0+Ql/Pl*Rtv 

52  =  R4+Ql/Pi*R<: 

PDN  *  -<Sl*R3+S2-*Pl  )/(Si+S2«R5) 

T  =  PDN*ALPHA  ( 2  ■'  ■'  ALPHA  ( 1 )  +  1.0 
S  -  PCN*EPS2/£P2t*ALFHA< 1 >/AlPHA<2>  +10 
PI  =  ALPHA ( 1  >+ALPHA(  1  ) 

V  -  SAMMA2+S/P1  -  T+S0EETA/EFS1/AK02 
P  -  OAHMA* < T/EPS! /AK02+S/P 1 */V 
01  =  <  GAMMA2+S03ETA)  /  ALPHA ( 1  •  /BETA! N )  /0I1EQA/V 
Q  *  01/EPS1/ZEP 
R  *  '31*S*T/ZMU 

11  =  Jl  +  1 

12  *  J1+J2 

13  =  J1+J2+1 

DO  oOO  1*1. J1 

600  SC  i  I !  »  COS ( BETA  • N ) *XP ( I > ) 

DO  610  1=11, 12 

610  SC<I>  =  SIN( BETA (  N ! *XP ( I ) > 

DO  620  1  =  13.  NP 

620  SC(I>  =  COS (BETA' N>*XP(I)> 

XSCJ  =  0.  O 
DO  50  J=l.  J1 

XSCJ1  =  SIN (BETA <N)+X(  J+l ) ) 

XT  «  XSCJ1  -  XSCJ 
XSCJ  *  XSCJ1 
DO  20  1=1, J1 

ZU.J)  «  2(1.  J)-G*XT*3C(I> 

20  CONTINUE 

DO  30  1  =  11,  12 

Z<I,J>  =  Z(I.  J)+P*XT*SC(I ) 

30  CONTINUE 

DO  40  1  =  13.  NP 

Z(I.J)  =  Z(I.  J)+0*XT*SC(I) 

40  CONTINUE 
50  CONTINUE 

XSCJ  =  COS (BETA>  N  )*X  ( !l)) 

DO  90  J=I1, 12 

XSCJ1  =  COS  <  BETA  <  N  >  *X ( J+l ) ) 

XT  *  XSCJ1  -  XSCJ 
XSCJ  =  XSCJ1 
DO  60  1=1,  J1 

Z(I.J)  =  Z(I.  J>-P*XT*SC(I) 

60  CONTINUE 

DO  70  1=11, 12 

Z(  I.  J>  =  Z<  I.  J)-R*XT+SC  ( I  ) 

70  CONTINUE 

DO  80  1  =  13.  NP 

Z(I.J)  =  Z(I. J)-F*XT+SC(I) 

80  CONTINUE 
90  CONTINUE 

XSCJ  =  SINCBETA  JJHM  13)  > 

DO  130  J=I3.NP 

XSCJ1  »  SIN(EETA'N>*X:.I+1)) 

XT  =  XSCJ1  -  XSCJ 
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XSCJ  =  X3CJ1 
DO  lOO  1  =  1.  J1 

Z(I  J>  =  Z< I, J>-*-u*XT*SC( 1 > 

100  CONTINUE 

DO  110  1=11, 12 

Ztl.J)  =  ZCI. J)-F*XT*SCi I) 

110  CONTINUE 

DO  120  1=13, NP 

Z(I,U>  =  ZCI,  J)-rO*XT*SC<  I  > 

120  CONTINUE 
130  CONTINUE 
200  CONTINUE 
C 

IE  (  I  IMP  .EG.  10  >  00  TO  700 

CALL  CCECO  (Z,  LDZ. NP.  IPVT, RCCND. WORK) 

IF  (  I  IMP  EG.  1  )  40  TO  400 

CALL  COED  I  ( Z . LDZ . NP ,  I PUT , DET , WORK,  101 
IF  <  IFACT  NE.  0  >  00  TO  32S 

FACTOR  «  R£AL<D£T;2> I 
IFACT  =  10 

325  DETERM  =  DET(l)  *  10. O^* (REAL. I DET ( 2 ) >  —FACTOR ) 

GAM  »  REAL  (  G4MM4 ) /AKO 
ATT  =  -AIMAG(GAMMA)+DBN 
EDET  =  REAL ( DET .  2 )  ) 

WRITE  (6,330)  GAM  ,  ATT  .  DET(1>  ,  EDET  ,  RCOND 
330  FORMAT  I.  1H  ,  4X.  F8.  4.  2X,  Ell  4,  3X.  F7.  4,  1  i,  F7.  4,  IX,  *E», 
1  FS  0.  2X,  E10  3  > 

RETURN 

C 

400  CALL  CGESL  <  Z,  LDZ,  NP.  IPVT.  WORK,  0) 

P  ■  (0.  0,  0.  0) 

Q  =  (O.  O.  0  0) 

DO  SCO  1=1, J1 

P  =  P— WORK  ( I )  *  ( >.  i  1  +  1 )  — X  ( I )  ) 

500  CONTINUE 

DO  510  1=11. 12 

a  =  a+woRK<i)*(<(i+n-xm> 

510  CONTINUE 

WRITE  < 1,450) 

450  FORMAT  <  ///  ) 

DO  520  1  =  1- NP 

WORK  I  I  )  =  WORK  (  I  ;•  /Q 

FLDR  =  REAL  <  WORK ( I ) ) 

FLDI  =  A I MAG  <  WORK  < I ) ) 

FLDA  =  CABS ( WORK ( I ) ) 

FLDC  =  PHASE  < WORK < I i > 

WRITE  (1,470)  /F ( I >. FLDR. FlDI,  FLDA,  FLDP 
470  FORMAT  (  5X,  'X  =  ' ,  E9  3,  3X.  '  (REAL)  ' ,  E13  7,  5X. 

1  '<  IMAGINARY)  ',£13  7,  5X,  '(ABSOLUTE)  El  3.  7. 

2  5X,  '(PHASE)  '.  F7.  2  > 

520  CONTINUE 

ZIMP  =  G/P/2.  0 
RETURN 
C 

700  LOU  =  1 
LDV  =  1 
JOB  =  0 
C 

CALL  CSVDC  (Z.  LD7,  NP- NP.  SC.  EA,  U.  LDU.  V.  LDV 
1  WQR».  JOB,  INFO; 

RCOND  =  REAL  <  SC  "JP  )  /SC '  1 )  ) 

GAM  =  REAL  uImMMA ) / AKO 
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ATT  =  -AIMAG(GAMMA)*CBN 

WRITE  <6.  710>  '-AM-  ATT,  RCOND.  INFO 

’10  FORMAT  (  1H  .  5X  GAM- ' ,  F8  4  OX,  'ATT-  ' ,  Eli;  n,  3X, 
1  ROUND-'. £10.  4,  3X.  'INFO-'.  12  > 


SUBROUTINE  POIN"a  ( X,  NP>  Jl.  J2>  U3,  Q) 

GENERATES  NON-UNIFORM  DISTRIBUTION  OF  THE 
MATCHING  POINTS. 

COMMON  /BL0C40/A,  B.  W,  Dl,  D2,  D3.  AKO.  PI,  ZMU,  2EP,  OMEGA 
DIMENSION  XUOl  ; 

YO  -  EXP(-G) 

Yl  -  EXP<Q+<A-B',A/2.  O) 

Y2  -  EXP<a*<B-U,,'4i 
YT  »  4.  O-YO-2.  0»Y  1-V2 

11  -  1 NT < FLOAT (Nr  Nil  O-YO  > ,  VT+O.  5  > 

12  -  INT < FLOAT <NF  >*( 1 .  0— Yl  > /  VT+O.  5> 

13  -  NP - I 1 -2* I 2 

X  <  1  >  -  0  0 
DO  10  1-2, 11 

X(I)  -  ALOGl FLOAT  <  1-1  )-*  ( 1.  O-YO)  /FLOAT  (Ill*  Y0)*A/Q  ♦  A 
10  CONTINUE 
X<I1+1)  -  A 
DO  20  1-2, 12 

XUl  +  I)  »  -ALQGi  1 .  O-FLOAT  ( I-i  )*<  1.  O-Yl )  /FLOAT  (12)  )*A/Q  +  A 
20  CONTINUE 

X  ( 1 1+12*1 )  -  (A-Sj/2.  0 
DO  20  1-2,  12 

XU1  +  I2+I)  -  ALC‘0  <  FLOAT  <1— 1)*(1.  0— Yl )  /FLOAT  ( 12) +Y1  ) 

1  *A/u  +  E 

30  CONTINUE 

X< 11+2*12+1 )  -  2 
DO  40  1-2, 13 

X< 11+2*12+1,  -  -4L0G< 1.  0-FL0AT< 1-1 )*( 1. 0-Y2) 

1  /FLOAT (13)) *A/0  +  B 

40  CONTINUE 
X  <  NP  :  1 )  =  W 
Jl  -  II 
J2  -  2*12 
J3  »  13 
RETURN 
END 


SUBROUTINE  PERIOD  i AKO. L, ZC ■ PR , GAM, SHIFT . ATTEN ) 

CALCULATES  COMPLEX  PROPAGATION  CONSTANTS 
OF  THE  PEP  IODIC  STRUCTURE. 
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COMPLEX 

REAL 


-liaatiiaMBVaial— 

-1-  -2- 
•*»a*aaa»»aa«a*— 

< - Ll - >  < - L2 - > 

X  ZC  (2  •  IK/ PR  <  2 ) i GAMMA, Z 

L<2).  Ll.  L2.  BETA <2 .  ATT <2) 


DEN  -  8  A859E-3 
Ll  ■*  L ( 1  > 
l2  *  L<2> 

ZK  *  ZC(1)/ZC<2- 

Z  »  ( <Zk+l.  0)**2-»CC05'>FR<  1  )-»t-l+PR»2>*L2.' 

1  -<ZK-1.  0>**2-*GCOS(PR<1>*L1-PR<2>*L2>  .<  •  ■*  O/ZK 

gamma  =  CLQG< Z+GSGRT <  Z*Z-1 > • >  (L1+L2  > 

SHIFT  »  6.  283185307179/ (L1+L2) /AKO 

IF  >  A I MAG (GAMMA'  .  LT.  0.0  )  GAMMA  «  -GAMMA 

GAM  =  AIMAG(GAMMA)/AKO 

ATTEN  *-  REAL ( GAMMA )*DBN 

RETURN 

END 


—  END  OF  THE  PROGRAM  — 
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